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This  scientific  report  presents  the  results  of  a  dedicated  experiment  that  was  conducted  within  the  framework  of  the  Low-latitude 
ionospheric  Sensor  Network  (LISN)  observatory  to  measure  the  characteristics  of  medium-scale  (hundreds  of  km)  Traveling 
Ionospheric  Disturbances  (TIDs)  as  they  transit  through  the  low-latitude  ionosphere.  A  small  array  of  3  GPS  receivers  separated 
by  4-5  km  placed  in  a  triangular  configuration  was  installed  near  Huancayo  in  Peru  possessing  several  characteristics  of  a  radio¬ 
interferometer.  During  the  campaign  days,  17-30  July  2008,  TIDs  were  observed  daily.  On  July  20,  2008  between  22  and  24  UT 
several  TIDs  moved  across  the  small  array  of  GPS  receivers  with  a  velocity  near  130  m/s,  were  directed  northward  and  had 
wavelengths  close  to  450  km.  Other  GPS  receivers  that  were  operating  hundreds  of  km  away  from  Huancayo  show  also  similar 
TEC  traces  and  provide  a  phase  velocity  equal  to  150  m/s.  This  value  was  measured  using  the  GPS  at  Piura,  Cuzco  and  Huancayo. 
Based  on  this  positive  result,  we  conclude  that  small  and/or  regional  arrays  of  GPS  receivers  can  be  used  at  low  latitudes  to  study 
the  role  that  gravity  waves  may  have  on  seeding  plasma  bubbles. 


1.  Introduction 

Traveling  ionospheric  disturbances  (TIDs)  are  the  iono¬ 
spheric  manifestation  of  neutral  density  oscillations  called 
atmospheric  gravity  waves  (AGWs  or  GWs).  AGWs  originate 
in  the  troposphere  or  stratosphere,  and  often  propagate 
to  F-region  heights,  where  they  induce  oscillations  in  the 
plasma  density  via  ion-neutral  coupling.  When  AGW  winds 
drive  plasma  along  magnetic  field  lines  via  ion-neutral 
collisions,  “classical”  TIDs  result  [1-3].  However,  at  night 
the  electrodynamics  are  favorable  to  the  production  of 
“nonclassical”  or  “electrodynamic”  TIDs,  such  that  vertical 
movement  of  the  F-layer  results  [4-6].  TIDs  that  reach 
F-region  altitudes  are  divided  into  two  classes:  medium- 
scale  (MSTIDs)  and  large-scale  (LSTIDs).  MSTIDs  primarily 
originate  from  tropospheric  weather  systems  and  mountain 
turbulence;  depending  on  background  wind  conditions  and 
other  factors,  these  may  propagate  in  any  direction.  LSTIDs 


are  produced  by  auroral  energy  injection  during  geomag¬ 
netic  storms.  They  typically  propagate  meridionally  from 
polar  regions  towards  the  equator.  AGWs  have  a  wide  range 
of  observed  wavelengths  and  propagation  velocities,  fre¬ 
quently  having  horizontal  wavelengths  of  up  to  hundreds  of 
kilometers  and  speeds  of  hundreds  of  meters  per  second.  As 
they  propagate  from  lower  altitudes,  their  amplitudes  grow 
exponentially.  At  mid-latitudes,  AGWs  have  been  observed  to 
have  preferred  propagation  azimuths  that  change  with  season 
[7];  it  is  logical  to  assume  that  similar  seasonal  patterns  exist 
at  low  latitudes. 

Gravity  waves  are  widely  believed  to  play  a  critical  role  in 
seeding  the  development  of  equatorial  plasma  depletions  in 
the  nighttime  low-latitude  ionosphere.  Current  understand¬ 
ing  holds  that  both  favorable  ionospheric  conditions  [8]  and 
an  initial  seed  perturbation  on  the  bottomside  are  necessary 
to  create  rising  depletion  channels.  If  both  are  present, 
the  Rayleigh-Taylor  instability  creates  intense  polarization 
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electric  fields  within  the  perturbations,  causing  the  depleted 
plasma  inside  to  E  X  B  drift  upward  into  the  topside  [9,  10]. 
Gravity  waves,  [11],  as  well  as  penetration  electric  fields, 
dynamics  associated  with  the  post-sunset  vortex  [12],  and 
E-region  structure  mapped  to  F-region  altitudes  along  mag¬ 
netic  held  lines  [13]  have  all  been  cited  as  possible  seeding 
mechanisms. 

To  unravel  the  importance  of  the  gravity  wave  seeding 
mechanism,  it  is  important  to  observe  and  characterize 
both  early-evening  gravity  waves  and  the  equatorial  plasma 
depletions  seen  later  in  the  evening.  However,  most  current 
observing  techniques  have  severe  limitations.  Ionosondes 
can  measure  the  GW  vertical  period,  but  they  are  not  sen¬ 
sitive  to  other  parameters.  A  steerable-beam  radar  such  as 
Altair  can  measure  the  large-scale  zonal  wavelength  and  am¬ 
plitude,  but  the  GW  must  be  traveling  roughly  orthogonal  to 
the  radar  beam  to  be  detected.  Phased-array  radar  systems 
such  as  the  Japanese  MU  system  [14,  15]  are  sensitive  to 
GW  horizontal  wavelength,  propagation  azimuth,  and  am¬ 
plitude,  but  these  systems  are  expensive  to  deploy  and 
operate.  01  (557.7  nm,  emission  altitude  ~96km)  and  OH 
band  (720-910 nm,  ~86km)  airglow  imagers  can  measure 
GW  horizontal  wavelengths,  propagation  azimuths,  and 
amplitudes,  but  at  altitudes  far  below  the  bottomside  [16, 
17].  OI  (630  nm  emission)  optical  measurements  have  been 
used  to  characterize  the  seasonal  and  solar  cycle  variability  of 
MSTIDs  [18,  19].  OI  has  an  emission  layer  between  roughly 
200-300  km,  which  is  very  close  to  the  bottomside  [20].  But 
even  here,  airglow  imagers  operate  only  during  the  nighttime 
and  require  good  seeing  conditions  (no  clouds)  to  operate. 

The  radio  beacon  technique  [7,  21-23]  overcomes  many 
limitations  of  other  observation  methods,  enabling  simulta¬ 
neous  characterization  of  key  GW  parameters.  Here,  satellite 
radio  beacon  signals  are  measured  by  stationary  ground 
receivers;  TIDs  are  evident  as  fluctuations  in  total  elec¬ 
tron  content  (TEC)  derived  from  the  beacon  signals  [24]. 
Although  TEC  is  an  integrated  quantity,  various  processes, 
including  wind-filtering  [7]  and  breaking  [25],  prevent 
AGWs  from  penetrating  far  into  the  thermosphere.  Obser¬ 
vations  [26]  and  model  results  [27]  have  shown  that  most 
TIDs  do  not  rise  above  300  km.  Thus,  the  bulk  of  the  TEC 
perturbations  occur  on  the  F-region  bottomside  and  the 
upper  E-region.  Moreover,  after  the  E-region  decays  in  the 
evening,  TEC  fluctuations  become  highly  localized  in  the 
bottomside  region.  Thus,  while  radio  beacon  TEC  measure¬ 
ments  are  integrated  over  the  whole  altitude  range  of  the 
ionosphere,  the  narrow  vertical  localization  of  GW-induced 
TEC  fluctuations  makes  this  technique  a  powerful  tool  for 
studying  early-evening  TIDs.  In  the  radio  beacon  technique, 
also  known  as  radio-interferometry,  phase  differences  mea¬ 
sured  at  the  various  stations  are  used  to  determine  TID 
velocity,  propagation  azimuth,  and  amplitude. 

Recently,  the  radio-interferometric  technique  has  been 
adapted  for  use  with  GPS  satellites  [28-30] .  This  new  innova¬ 
tion  makes  it  possible  to  utilize  inexpensive,  easily  deployed 
GPS  receivers  to  study  gravity  waves  at  a  wide  variety  of 
locations.  The  large  number  of  GPS  satellites  in  orbit  makes  it 
possible  to  continuously  monitor  the  bottomside  region  over 
a  given  receiver  array. 


Figure  1:  Geographic  location  of  the  Huancayo  radio  interferome¬ 
ter  test  array. 


This  paper  describes  measurements  observed  with  a 
three-station  array  of  GPS  receivers  that  were  located  near 
Huancayo,  Peru.  These  receivers  constitute  the  first  low- 
latitude  radio-interferometer  array  dedicated  to  studying 
TIDs.  The  hardware  layout  and  the  analysis  technique,  that 
was  used  to  derive  the  TEC  perturbation  associated  with 
the  TIDs,  are  described  in  the  next  section.  The  phase 
velocities,  scale  sizes,  and  the  direction  of  motion  of  the  TIDs 
calculated  using  the  SADM-GPS  methodology  are  presented 
in  Section  3.  The  TIDs  wave  characteristics  derived  using  a 
cross- correlation  method  are  introduced  in  Section  4.  The 
cross- correlation  method  was  also  applied  to  TEC  mea¬ 
surements  recorded  by  other  low-latitude  ionosphere  sensor 
network  (LISN,  http://lisn.igp.gob.pe/)  GPS  receivers  that  are 
operating  hundreds  of  kms  apart  from  Huancayo.  The  results 
of  this  analysis  are  described  in  Section  5.  Discussion  and 
conclusions  follow  in  Sections  6  and  7. 

2.  Observations 

Since  the  late  90s  several  publications  have  reported  the 
characteristics  of  mid-latitude  GWs  using  the  GPS  radio¬ 
interferometry  (GPS-RI)  method  [28,  30].  Based  on  this 
success,  we  conducted  a  campaign  in  South  America  aiming 
to  characterize  the  local  time,  duration,  and  wave  properties 
of  GWs  that  are  transiting  near  the  magnetic  equator  in 
the  Peruvian  sector.  We  installed  two  GPS  receivers  near  the 
city  of  Huancayo  (12.042°S;  75.321°W)  due  to  its  proximity 
to  the  magnetic  equator,  the  existence  of  a  permanent  GPS 
receiver  in  that  city  and  the  presence  of  nearby  towns.  The 
additional  GPS  receivers  were  temporarily  installed,  one  at 
the  town  of  Chupaca  (12.061°S;  75.292°W)  and  the  other  at 
Sicaya  (12.021°S;  75.282°W).  All  3  GPS  receivers  of  the  small 
network  of  GPS  receivers  consist  of  GSV4004B  Novatel  GPS 
using  geodetic  antennas.  Similarly,  Novatel  GPS  receivers 
were  also  installed  in  the  cities  of  Piura  and  Cuzco.  However, 
the  receiver  in  Iquitos  consists  of  a  CRS1000  LEICA  GPS  that 
uses  a  choke  ring  antenna.  Figure  1  indicates  the  geographic 
locations  of  the  three  sites,  where  GPS  receivers  operated 
between  July  17 — 30,  2008  and  July  30.  As  indicated  in 
Figure  1,  the  distance  between  GPS  receivers  was  between 
3.9  and  4.8  km.  They  were  distributed  forming  an  isosceles 
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Figure  2:  TEC  versus  UT  measured  by  the  Huancayo  GPS  receivers;  red,  green,  and  blue  traces  represent  TEC  values  from  the  three  receivers; 
black  lines  represent  the  estimated  background/unperturbed  TEC  values.  Each  of  the  four  panels  represents  signals  from  one  of  four  GPS 
satellites  (PRN  =  07,  14,  22,  and  25). 


triangle  to  provide  good  estimates  of  the  phase  velocity  as  a 
function  of  both  latitude  and  longitude.  Due  to  the  relatively 
short  distance  between  receivers  and  the  relatively  large  phase 
velocity  of  GWs  that  can  vary  between  100  and  500  m/s  [31], 
we  opted  to  sample  the  phase  and  group  delay  of  both  LI  and 
L2  frequencies  once  per  second  ( 1  Hz). 

Each  of  the  four  frames  of  Figure  2  present  TEC  values 
measured  by  all  3  GPS  receivers  that  operated  near  Huancayo 
on  July  20,  2008.  The  TEC  values  from  different  receivers 
have  been  color  coded  to  differentiate  them  and  to  indicate 
their  similarity.  The  TEC  values  were  computed  following 
the  standard  procedure  that  involves  using  the  differential 
phase  and  pseudorange  values  of  both  LI  and  L2  frequency 
signals  [32]  and  leveling  the  differential  phase  to  the  average 
differential  code  values  to  eliminate  phase  ambiguity.  The 
absolute  values  of  TEC  were  then  calculated  by  introducing 
the  differential  satellite  biases  published  by  the  University 
of  Bern  and  calculating  the  receiver  bias  by  a  minimization 
process  of  the  TEC  variability  between  0200  and  0600  local 
time.  The  thin  black  line,  displayed  in  each  of  the  4  frames 
of  Figure  2,  indicates  an  estimate  of  the  background  or  un¬ 
perturbed  TEC  values.  This  estimate  takes  into  account  the 
typical  daily  variability  of  the  TEC  values  that  is  a  function 
of  the  local  time  variability  produced  by  solar  EUV  radiation, 
changes  in  tides,  presence  of  external  electric  fields,  and  more 
importantly  the  TEC  variability  introduced  by  the  motion 
of  the  GPS  satellites  that  continuously  samples  different  lati¬ 
tudes  and  longitudes.  These  background  TEC  values  are  used 


to  estimate  the  TEC  perturbation,  or  in  other  words,  derive 
the  contribution  produced  by  the  passage  of  the  gravity 
waves.  The  TEC  curves  corresponding  to  each  satellite  pass 
were  fitted  to  a  4th-order-polynomial  to  filter  out  the  diurnal 
variability  and  to  remove  spatial  changes  in  TEC  due  to 
the  GPS  satellite  motion  and  to  evaluate  TEC  perturbations 
(TECP)  containing  time  scales  of  3  hours  or  less.  Similar 
analysis  was  conducted  by  Valladares  et  al.  [24]  to  estimate 
the  TEC  perturbation  attributed  to  the  passage  of  large-scale 
gravity  waves.  Equatorial  plasma  bubbles  can  also  produce 
significant  perturbations  in  the  TEC  traces.  However,  their 
effect  consists  of  rapid  and  deep  negative  excursions  of  the 
TEC  traces,  commonly  called  TEC  depletions.  To  avoid  any 
possible  erroneous  identification,  the  GW/TID  campaign 
was  carried  out  during  a  period  of  minimum  bubble  activity, 
in  which  no  TEC  depletions  were  observed.  As  shown  in  the 
following  figures,  all  the  analysis  methods  are  based  on  the 
difference  between  the  measured  and  the  estimated  back¬ 
ground  “unperturbed”  TEC  values.  This  differential  value  is 
called  TEC  perturbations  (TECP). 

Figure  3  shows  TECP  values  measured  on  3  different 
days.  The  top  panel,  corresponding  to  July  17,  2008,  displays 
perturbations  that  reach  1  TEC  unit  between  20  and  22  UT. 
Smaller  amplitude  TEC  perturbations  (±0.4  TEC  units)  are 
observed  between  13  and  19  UT.  The  noise  level  in  this  and 
other  panels  is  about  ±0.2  TEC  units.  The  TEC  perturbations 
of  the  center  panel  (July  20, 2008)  exceed  1  TEC  unit  between 
22  and  24  UT  and  exhibit  moderate  levels  (±0.6  TEC  units) 
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Figure  3:  TECP  versus  UT  at  Huancayo  for  three  days  in  July.  Data  from  various  GPS  satellite  signals  are  plotted;  each  satellite  is  distinguished 
by  a  different  color. 


of  TEC  perturbations  between  14  and  16  UT.  Note  the 
different  amplitude  scale  that  is  used  in  this  panel.  The 
bottom  panel  corresponding  to  July  15,  2008  depicts  two 
large  patches  with  significant  levels  of  TEC  perturbations. 
The  first  occurs  during  early  evening  hours  (00-04  UT)  and 


the  second  in  the  afternoon  hours  (19-24 UT).  In  sum¬ 
mary,  during  the  campaign  TEC  perturbations  are  observed 
preferentially  a  few  hours  before  sunset  (24  UT).  A  second 
interval  of  TEC  perturbations  occurs  at  different  times  of  the 
day. 
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3.  Analysis  of  TID  Characteristics  Using 
SADM-GPS  Algorithm 


We  have  calculated  the  TIDs  horizontal  phase  velocity  14(f) 
and  the  azimuthal  propagation  direction  of  the  phase  front 
a{t)  using  the  statistical  angle  of  arrival  and  doppler  method 
for  GPS  interferometry  (SADM-GPS)  developed  by  Afraim- 
ovich  et  al.  [28-30].  This  method  assumes  that  the  iono¬ 
sphere  changes  in  a  nondispersive  manner  due  to  the  passage 
of  the  TID  or  a  packet  of  TIDs.  Under  this  assumption, 
Afraimovich  et  al.  [28,  30]  calculated  the  velocity  and  direc¬ 
tion  of  propagation  of  mid-latitude  TID  packets  during  quiet 
ionospheric  conditions.  The  SADM-GPS  represents  the  TEC 
perturbations  by  the  form:  I{x,y,t )  =  d(f)sin(Of  -  kxx  - 
kyy  +  (po).  This  expression  implies  that  a  monochromatic 
or  a  largely  dominant  TID  propagates  through  the  field-of- 
view  of  the  small  receiver  network.  For  events  when  more 
than  one  TID  is  propagating  through  the  measured  volume, 
this  representation  will  be  invalid,  and  a  spectral  characteri¬ 
zation  of  the  TIDs  is  more  appropriate.  Using  the  equations 
proposed  by  Afraimovich  et  al.  [28],  it  can  be  demonstrated 
that  the  TID  wave  parameters  can  be  calculated  using: 


a(t)  =  arctan(|g) 


j  =  arctan(lj, 

(1) 

4'(f) 

m ’ 

(2) 

m 

ry{ty 

(3) 

u(t )  = 


UX(t)Uy(t) 


1/2  ’ 


(4) 


Vk{t)  =  u(t)  +  wx{t )  sin  a(t)  +  wy(t )  cos  a(t), 


(5) 


Gx 


(Ta(TECP£  -  TECPc)  -  Yc( TECPB  -  TECPA)) 
XAYC-XCYA 


(6) 


r  _  (Xc(TECPb  -  TECPa)  -  XA(TECPa  -  TECPc)) 

'  xAyc-xcyA 

(7) 

where  I'x(t),  I' y{t),  and  /'t(t)  are  the  spatial  and  time  deriv¬ 
atives  of  the  TECP  values.  ux  ( t)  and  uy  ( t )  are  the  propagation 
velocities  of  the  phase  front  along  the  x-(east)  and  y-(north) 
axes  in  a  frame  of  reference  centered  at  the  Huancayo  station. 
wx(t)  and  wy(t)  are  the  x  and  y  projections  of  the  sub-iono¬ 
spheric  intersection  point  velocity.  Note  that  Gx  and  Gy  sym¬ 
bols  used  in  Figure  4  are  equivalent  to  the  I'x(t)  and  I' y(t) 
derivatives.  As  the  velocity  of  the  subionospheric  intersection 


(80-100  m/s)  is  comparable  to  the  speed  of  the  GWs,  it 
becomes  imperative  to  subtract  this  effect  from  the  estimated 
phase  velocity  u{t).  Equation  (5)  includes  this  correction  in 
which  the  component  of  the  satellite  velocity  in  the  direction 
perpendicular  to  the  phase  front  is  subtracted  to  determine 
the  true  horizontal  velocity  of  the  TID.  The  symbols  XA,  YA, 
Xc,  and  Y(  are  defined  in  Figure  1  and  correspond  to  the 
coordinate  distances  of  Sicaya  (A)  and  Chupaca  (C)  in  a 
Cartesian  system.  To  eliminate  noise  fluctuations  containing 
periods  less  than  5  min,  a  low-pass  filter  was  applied  to  all 
TECP  traces  shown  in  subsequent  figures. 

Figure  4  presents  the  results  of  applying  (1)  through  (7) 
to  the  TECP  values  measured  on  July  20,  2008  between  21 
and  24  UT.  This  figure  follows  closely  the  format  presented 
by  Afraimovich  et  al.  [28]  in  Figure  9.  Panels  (a)  and  (e)  show 
the  time  derivative  ( Gt )  of  the  TECP  measured  at  Huancayo 
from  GPS  satellites  13  and  25,  respectively.  Panels  (b)  and  (f) 
display  the  east-west  (Gx)  spatial  derivative  using  a  solid  line, 
and  the  north-south  ( Gy )  spatial  derivatives  using  dots.  Note 
that  (6)  {Gx)  and  (7)  ( Gy )  are  a  function  of  the  TECP  values 
measured  by  all  3  GPS  receivers  belonging  to  the  small  array. 
Panels  (c)  and  (g)  present  the  angle  of  the  wave  front,  a(f), 
that  was  calculated  using  (1).  It  is  noted  that  this  angle  is 
measured  clockwise  with  respect  to  geographic  north.  Panels 
(d)  and  (h)  introduce  the  phase  velocity  of  the  TID  calcu¬ 
lated  using  (5).  The  large  scatter  of  the  a(t)  and  14(f)  values 
is  caused  by  the  nonmonochromatic  nature  of  the  TECP  val¬ 
ues.  However,  the  average  lines  in  these  frames  indicate  our 
best  estimate  of  the  TID  phase  velocity  that  varies  between 
3°  and  6°.  The  estimation  of  the  horizontal  velocity  of  the 
dominant  TID  falls  between  107  and  143  m/s.  Comparison  of 
the  Gt  traces  for  satellites  13  and  25  (panels  (a)  and  (e))  show 
largely  different  magnitudes;  this  is  due  to  the  different  look 
angle  direction  that  varies  between  these  two  GPS  satellites. 
GPS  satellite  13  passed  200  km  away  from  Huancayo  and 
GPS  25  went  over  Huancayo  almost  overhead.  Mercier  and 
Jacobson  [33]  concluded  that  the  largest  TECP  values  were 
found  when  the  line  of  sight  made  a  small  angle  with  the 
phase  front.  The  satellite  trajectory  of  GPS  25  was  almost 
parallel  to  the  propagation  direction  of  the  phase  front. 

Figure  5  shows  hourly  histograms  of  the  statistics  of  the 
phase  front  angle,  a(t),  and  the  horizontal  velocity,  14(f), 
comprising  all  GPS  satellites  that  passed  near  Huancayo 
between  21  and  24  UT.  A  total  of  6  GPS  satellites  observed 
significant  TEC  perturbations  associated  with  the  TID.  Each 
frame  of  Figure  5  contains  a  vertical  arrow  to  point  out  the 
location  of  the  statistical  average.  The  propagation  angle 
of  the  TIDs  (left  frames)  indicates  a  consistent  northward 
direction.  The  phase  velocity  of  the  TIDs  varies  between  60 
and  120  m/s  with  the  maximum  value  occurring  between  22 
and  23  UT. 


4.  Analysis  of  TID  Characteristics  Using 
Cross-Correlation  Algorithm 

This  section  describes  the  cross-correlation  method  (CCM) 
and  presents  our  results  of  the  TID  phase  velocity,  its 
direction  of  propagation  and  scale  size.  The  aim  of  the  CCM 
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Figure  4:  Derivation  of  different  parameters  used  by  the  SADM-GPS  method  for  GPS  satellites  13  and  25.  Panels  (a)  and  (e)  show  Gt;  (b) 
and  (f)  Gx  solid  curve  and  GXl  dots;  panels  (c)  and  (g)  show  the  direction  of  propagation  a(t);  (d)  and  (h)  present  the  horizontal  velocity, 
V*(f). 
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Figure  5:  Distributions  of  the  azimuth  angle  on  the  left  and  velocities  on  the  right,  corresponding  to  the  calculated  dependencies  of  a(t)  and 
14(f).  The  large  arrows  indicate  the  location  of  the  mean  value  for  each  histogram. 


method  is  to  calculate  the  running  cross-correlation  function 
using  TECP  traces  originated  from  any  of  the  GPS  satellites, 
but  measured  by  two  different  receivers.  The  cross-correla¬ 
tion  function  is  calculated  every  minute  using  segments  of 
TECP  data  one  hour  long.  We  use  TECP  values  from  the 
same  GPS  satellite  due  to  the  motion  of  the  satellite  that 
produces  the  same  temporal  delay  or  advance  on  receivers 
placed  at  relatively  close  distances.  This  satellite  motion  effect 
can  be  removed  after  we  calculate  the  phase  velocity.  If  signals 
from  different  GPS  satellites  are  used,  it  would  be  necessary 
to  remove  the  effect  of  the  satellite  motion  before  the  CCM 
method  is  applied.  Figure  6  shows  the  results  of  the  CCM 
method  based  on  signals  from  the  GPS  13  satellite  recorded  at 
Huancayo  and  Chupaca.  The  right  frame  of  Figure  6  shows 
the  same  analysis  using  the  TECP  values  measured  at  Huan¬ 
cayo  and  Sicaya.  The  lower  panels  show  the  TECP  traces  for 
GPS  13  in  TEC  units  observed  on  July  20,  2008  and  between 
20  and  24  UT.  The  lower  left  frame  displays  the  TECP  values 
measured  at  Huancayo  (blue  line)  and  Chupaca  (green  line). 
Due  to  the  coarse  resolution  of  the  plot,  the  small  time  delay 
between  the  traces  becomes  invisible.  The  cross-correlation 
functions,  displayed  in  the  upper  frames,  indicate  that  the 
time  delay  is  of  order  10  sec.  Note  the  color  scale  of  the  cross¬ 
correlation  functions  varies  between  0.9  and  1.0  due  to  the 
high  degree  of  coherence  of  both  signals.  To  obtain  the  scale 


size  of  the  TECP,  we  fitted  an  ensemble  of  monochromatic 
sinusoids  to  each  segment  one-hour  long.  The  numbers 
printed  in  the  lower  frame  point  out  the  scale  size  in  minutes 
of  the  best  fit  for  each  hour  of  TECP  data.  These  results  are 
affected  by  the  motion  of  the  GPS  satellites  and  need  to  be 
corrected  by  subtracting  the  apparent  motion  of  the  sub- 
ionospheric  intersection  in  the  direction  perpendicular  to  the 
phase  front.  The  propagation  direction  and  velocity  of  the 
TID  was  obtained  using: 


a{t)  =  arctan 


YcTb-a  ~  YaTb-c 
XcTb-a  ~  XaTB-c 


Vh(t) 


{Yc  cos(a(f))  -  Xc  sin («(())) 
Tb-c 


+  wx(t)  sina(f)  +  wy{t)  cos  a(t). 


(8) 

(9) 


Tb-A  and  Tb-c  are  the  time  delay  between  the  Huancayo 
and  Sicaya  TEC  perturbations  and  the  offset  time  between 
the  Huancayo  and  Chupaca  TECP  traces.  As  in  the  previous 
section,  wx(t)  and  wy(t)  are  the  projections  of  the  sub-iono¬ 
spheric  intersection  point  velocity.  A  negative  offset,  as  seen 
on  the  upper  left  frame  indicates  that  the  TECP  observed 
at  Huancayo  lags  Chupaca’s  TECP.  A  positive  delay,  as  seen 
in  the  upper  right  frame  is  an  indication  that  Huancayo’s 
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Figure  6:  Cross-correlation  functions  (top  panels)  and  TECP  curves  (lower  frames)  for  GPS  satellite  13.  The  left  frames  show  the  CCFs 
for  TECP  signals  from  Huancayo  and  Chupaca.  The  right  panels  correspond  to  Huancayo  and  Sicaya.  The  cross-correlation  functions  are 
plotted  using  a  color  scale  that  varies  between  0.9  and  1.0.  The  black  line  that  crosses  near  the  middle  corresponds  to  the  maximum  and  it  is 
taken  as  the  time  delay  ( Tb-a  or  Tb-c)- 
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Figure  7:  Phase  velocity  of  the  TIDs  calculated  using  the  CCM 
method  and  based  on  the  data  collected  by  the  small  array  that  oper¬ 
ated  near  Huancayo. 


TECP  leads  Sicaya’s.  Based  on  this  timing  relationship,  it 
can  be  concluded  that  the  perturbation  associated  with  the 
TID  was  seen  first  at  Chupaca,  located  further  south,  then 
at  Huancayo  placed  near  the  center  of  the  array,  and  later 
at  Sicaya  placed  further  north  (Figure  1).  This  sequence  of 
events  indicates  that  the  TID  motion  was  in  the  northward 
direction. 

A  more  quantitative  estimation  of  the  TID  phase  velocity 
was  obtained  using  (8)  and  (9).  Figure  7  shows  the  magni¬ 
tude  and  direction  of  the  TID  phase  velocity  represented  by 
arrows.  The  size  of  the  arrow  is  proportional  to  the  phase 
velocity  and  the  azimuthal  rotation  of  the  vector  is  a  function 


of  the  propagation  angle.  An  arrow  pointing  to  the  top 
indicates  a  propagation  direction  toward  geographic  north.  It 
is  evident  that  between  22  and  24  UT,  the  TID  moved  toward 
geographic  north.  A  wavelength  equal  to  449  km  is  obtained 
using  the  true  phase  velocity  of  the  TID  and  the  velocity  of 
the  satellite  sub-ionospheric  intersection  point.  The  phase 
velocity  varied  between  84  m/s  and  144  m/s  near  24  UT. 
These  values  are  consistent  with  the  results  of  the  SADM- 
GPS  method  (Section  3). 

5.  Analysis  of  the  TID  Characteristics  Using  GPS 
Receivers  Spaced  Hundreds  of  Kilometers 

We  also  report  our  investigation  of  the  TID  morphological 
properties  using  GPS  receivers  that  belong  to  the  low-latitude 
ionospheric  sensor  network  (LISN)  that  operate  at  different 
cities  separated  by  hundreds  of  kms.  LISN  is  a  distributed 
observatory  designed  to  nowcast  the  state  and  dynamics  of 
the  low-latitude  ionosphere  over  South  America.  The  LISN 
observatory  is  comprised  of  GPS  receivers,  flux-gate  mag¬ 
netometers  and  vertical  incidence-pulsed  ionospheric  radar 
(VIPIR)  ionosondes.  A  unique  feature  of  the  observatory  is 
the  ability  of  all  the  instruments  to  relay  their  observables 
in  a  near  real-time  manner  via  Internet.  At  the  time  of  this 
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Figure  8:  Same  as  Figure  3,  but  for  three  of  the  LISN  stations.  From  top  to  bottom:  Piura,  Iquitos,  and  Cuzco.  These  TECP  values  were  ob¬ 
tained  on  July  20,  2008. 


paper  (November  2011),  44  dual-frequency  GPS  receivers 
were  part  of  the  LISN  network.  During  the  time  of  the  cam¬ 
paign  a  total  of  126  GPS  receivers — including  the  LISN  GPS 
receivers — operated  in  the  South  American  continent.  TEC 
perturbations  from  these  receivers,  especially  the  ones  lo¬ 
cated  close  to  Huancayo,  were  used  to  provide  supportive 


evidence  for  the  TID  dynamic  characteristics  and  give  a  re¬ 
gional  view  of  the  presence  of  TIDs  over  South  America. 

Figure  8  shows  the  TECP  values  measured  on  July  20, 
2008  at  Piura  (5.17°S;  80.63°W),  Iquitos  (3.77°S;  73.27°W) 
and  Cuzco  (13.52°S,  71.96°W).  The  TEC  perturbations  in 
Piura  were  seen  mainly  between  21  and  24  UT.  Iquitos, 
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Figure  9:  Same  as  Figure  6,  but  for  CCF  between  Huancayo  and  Cuzco  on  the  left  side  and  Huancayo  and  Piura  on  the  right  side. 
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Figure  10:  Same  as  Figure  7,  but  for  TID’s  phase  velocities  that  were 
calculated  using  TECP  values  from  Huancayo,  Cuzco,  and  Piura. 


located  in  the  north-eastern  side  of  Peru,  registered  large 
TECP  values  between  16  and  24  UT.  Cuzco,  located  close  to 
the  magnetic  equator  and  250  km  east  of  Huancayo,  observed 
large  TEC  values  between  14  and  24  UT.  These  3  stations 
show  TECP  values  much  larger  than  the  noise  level  almost 
simultaneously.  However,  the  TECP  amplitude  is  different 
at  different  stations.  Figure  8  suggests  that  between  21  and 
24  UT,  TIDs  existed  in  the  western  region  of  South  America. 
Below,  we  discuss  further  this  Ending. 

The  lower  panels  of  Figure  9  display  the  TECP  values  for 
Huancayo  (in  blue),  Cuzco  (left  panel  in  red),  and  Piura 
(right  panel  in  red).  The  upper  left  frame  presents  the  CCFs 
calculated  using  the  Huancayo  and  Cuzco  TECPs,  showing  a 
relatively  short  delay  that  varies  between  +10  and  -15  min. 


Note  that  due  to  the  large  distances  between  the  stations,  the 
time  delays  between  these  station’s  TECP  values  are  of  order 
of  tens  of  minutes.  The  right  upper  frame  of  Figure  9  displays 
CCFs  for  Huancayo  and  Piura.  The  time  delay  between  the 
TECP  values  for  these  two  stations  is  ~45min.  This  value 
agrees  quite  well  with  the  time  offset  that  is  observed  in  the 
lower  panel.  Visual  inspection  of  the  minimum  value  of  the 
blue  curve  (Huancayo)  observed  at  -2230  UT  and  the  mini¬ 
mum  value  of  the  red  trace  seen  at  -2315  UT  indicates  that 
the  Huancayo  TECP  values  lead  the  Piura  TECP  values  by 
45  min.  Note  that  CCFs  in  the  upper  panels  also  show  sec¬ 
ondary  maximum  with  much  smaller  amplitudes  (-0.6) 
and  spaced  roughly  50-55  min  away  from  the  primary  maxi¬ 
mum.  In  the  case  of  a  perfect  sinusoidal  oscillation,  this  spac¬ 
ing  would  be  equal  to  the  wavelength  of  the  sinusoid.  How¬ 
ever,  since  MS-TIDs  are  more  properly  viewed  as  traveling 
wave  packets  [30],  this  spacing  corresponds  to  the  scale 
size  of  the  TID  packet.  It  is  also  pointed  out  that  the  CCF 
amplitude  scale  was  limited  to  values  above  0.4  and  did  not 
include  negative  values  as  they  would  indicate  anticorrela¬ 
tions.  If  we  had  used  a  CCF  amplitude  scale  varying  between 
- 1  and  1 ,  a  negative  peak  would  have  been  observed  between 
the  primary  and  secondary  maxima. 

Figure  10  shows  the  TIDs  phase  velocity  that  was  ob¬ 
tained  using  (8)  and  (9)  and  the  time  delays  presented  in 
Figure  9.  The  similarity  of  the  angle  a(t)  that  was  derived 
based  on  measurements  by  GPS  receivers  separated  by 
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Figure  11:  TEC  perturbations  recorded  by  several  GPS  sites  that  operated  in  South  America  on  July  20,  2008.  The  red  traces  correspond  to 
the  TECP  values  after  subtracting  the  daily  variability. 


hundreds  of  kms  and  the  a(t)  calculations  based  on  GPS 
receivers  that  were  installed  only  few  kms  apart  is  striking. 
Figures  7  and  10  show  a(t)  to  be  directed  very  closely  in  the 
geographic  northward  direction.  However,  the  amplitude  of 
the  horizontal  phase  velocity  seen  at  2230 UT,  varies 

between  lOOm/s  in  Figure  7  and  160  m/s  in  Figure  10.  The 
difference  in  the  Vj,(f)  values  is  reduced  at  2300  UT,  when 
135  m/s  is  shown  in  Figure  7  and  147  m/s  in  Figure  10. 

The  TEC  perturbations  for  all  126  GPS  receivers  that 
operated  in  South  America  on  July  20,  2008  are  presented  in 
Figure  11.  These  GPS  receivers  provide  a  good  coverage  over 
the  whole  South  American  continent,  except  for  latitudes 
southward  of  40° S.  Figure  11  shows  TECP  values  in  red, 
centered  along  the  satellite  trajectory  (thin  black  line),  that 
have  amplitudes  higher  than  0.4  TEC  units.  This  baseline 
value  is  about  4  times  the  noise  level  of  TEC  fluctuations. 
The  TECP  traces  of  Figure  1 1  demonstrate  that  TIDs  were 
restricted  to  the  north-western  side  of  South  America  be¬ 
tween  2200  and  2300  UT.  Similar  Figures  corresponding  to 
earlier  times,  not  shown  here,  display  the  TECP  traces  to 
be  located  further  south.  They  indicate  that  the  TIDs  were 
propagating  northward  and  were  probably  originated  in  the 
northern  part  of  Chile  (lat.  =  20°S). 

6.  Discussion 

We  conducted  a  campaign  at  Huancayo  in  Peru  to  assess  the 
effectiveness  of  GPS  receivers  in  detecting  and  measuring  the 


morphology  and  characteristics  of  TIDs  that  circulate  at  low 
latitudes.  To  accomplish  this,  we  temporarily  installed  two 
additional  GPS  receivers  near  the  Huancayo  station  to  form 
3  closely  spaced  GPS  receivers  having  a  geometry  that  can 
be  used  as  a  GPS  radio-interferometer.  The  TIDs  that  were 
observed  during  the  campaign  can  be  classified  as  medium- 
scale  disturbances  (MSTIDs).  They  have  horizontal  phase 
velocities  between  100  and  250  m/s,  wavelengths  of  several 
hundred  kilometers  and  periods  between  15  and  60  min  [34], 
We  also  indicate  that  the  TIDs  observed  on  July  20,  2008  are 
associated  with  AGW,  that  are  likely  generated  in  the  lower 
atmosphere. 

To  derive  the  TID  characteristics  we  applied  two  different 
methods  (SADM-GPS  and  CCM)  and  used  two  different 
datasets.  The  SADM-GPS  is  a  statistical  method  that  pro¬ 
vides  the  correct  phase  velocity  and  the  angle  of  propagation 
of  the  TIDs.  However,  it  requires  averaging  the  derived 
parameters  for  one  hour  and  for  all  the  GPS  satellites  with 
look  angles  above  35°  elevation.  The  CCM  method  seems 
more  robust,  and  it  is  less  affected  by  the  fact  that  more  than 
one  AGW  may  be  circulating  across  the  volume  probed  by 
the  GPS  receivers.  A  more  detailed  comparison  of  the  results 
provided  by  the  SADM-GPS  and  the  CCM  methods  shows  an 
excellent  agreement  in  the  angle  of  propagation,  as  both  indi¬ 
cate  the  angle  of  propagation  to  be  equal  to  +3°  ±  5°.  Cal¬ 
culation  of  the  phase  velocity  with  the  SADM-GPS  method 
indicates  an  average  velocity  of  120  m/s  between  22  and  23 
UT.  The  CCM  method  is  more  precise  and  indicates  that 


12 


International  Journal  of  Geophysics 


the  phase  velocity  varies  between  84  and  130  m/s  during  the 
same  interval. 

We  have  also  applied  the  CCM  method  to  TECP  values 
collected  by  other  GPS  receivers  that  belong  to  the  LISN 
network  and  are  operating  at  Piura,  Cuzco,  and  Iquitos,  all 
of  them  in  Peru.  We  found  very  good  correlation,  in  excess 
of  0.9,  between  the  TECP  traces  corresponding  to  stations 
placed  hundreds  of  kms  away.  We  were  able  to  calculate  that 
the  angle  of  propagation  a{t)  was  close  to  0°  (northward 
propagation)  and  the  horizontal  phase  velocity  V/,(t)  was 
between  130  and  160  m/s.  The  latter  values  are  larger  than 
the  numbers  obtained  using  the  small  interferometer  array  of 
3  GPS  receivers.  However,  several  geometrical  considerations 
need  to  be  accounted  with  more  precision.  One  of  them  is 
the  fact  that  the  projection  of  the  apparent  satellite  motion 
over  the  propagation  direction  of  the  TID  may  vary  during 
the  45  min  travel  time  between  the  LISN  stations.  It  is  also 
possible  that  the  TID  may  not  be  a  plane  wave,  and  the 
wave  front  may  have  a  circular  shape.  Figure  1 1  shows  that 
the  region  populated  by  TIDs  seems  to  extend  over  a  large 
region.  TECP  images  obtained  prior  to  Figure  11  indicate 
that  the  TIDs  originated  near  latitude  =20°  S.  We  do  not 
know  if  they  were  produced  by  deep  convection  cells  or  were 
a  response  to  orographic  features  in  northern  Chile.  We 
believe  that  the  success  in  getting  good  agreements  is  partially 
due  to  the  scale  size  of  the  TIDs  that  was  over  400  kms.  This 
distance  is  larger  than  the  separation  between  GPS  sites. 
Further  studies  will  determine  the  preferred  scale  sizes,  the 
phase  velocity,  and  the  direction  of  propagation  of  the  TIDs 
using  all  the  GPS  receivers  that  are  operating  in  South 
America. 

The  GPS  radio-interferometry  (GPS-RI)  method  is  one 
of  the  most  effective  techniques  with  which  to  observe  GW/ 
TID  activity  at  low  latitudes.  TIDs  are  notoriously  difficult  to 
characterize  due  to  their  large  size.  However,  GPS-RI  enables 
characterization  of  their  large-scale  features  with  ease,  in¬ 
cluding  their  horizontal  wavelength,  phase  speed,  and  pro¬ 
pagation  azimuths.  Since  GPS-RI  is  based  on  radio-frequen¬ 
cy  measurements,  the  technique  is  not  limited  by  weather 
considerations;  data  may  be  taken  at  all  times,  unlike  with 
optical  observations.  GPS  receivers  are  much  less  expensive 
to  deploy  than  other  instruments  of  similar  or  lesser  sen¬ 
sitivity  to  TIDs,  including  phased-array  radars,  steerable- 
beam  radar  systems,  and  ionosondes.  Moreover,  an  extensive 
global  array  of  GPS  receivers  already  exists,  allowing  for 
widespread,  continuous  monitoring  of  TID  activity.  Never 
before  has  this  been  possible. 

Additional  observations,  when  combined  with  GPS-RI 
observations,  yield  an  even  more  complete  TID  characteriza¬ 
tion.  In  particular,  measurements  of  the  thermospheric  neu¬ 
tral  wind  may  be  used  to  calculate  TID  intrinsic  frequency 
by  the  relation:  a>  =  KniVhit),  —  Un),  where  co  is  the  intrinsic 
frequency,  Kh  is  the  TID  horizontal  wavenumber  of  the  TID, 
and  Un  is  the  neutral  mean  wind  in  the  direction  of  TID 
propagation.  Further,  if  the  intrinsic  frequency  is  known, 
TID  vertical  wavelength  may  be  calculated  by  using  one 
of  the  dispersion  relations  [35,  36].  Thus,  a  very  complete 
picture  of  TIDs  results  with  only  GPS-RI  and  thermospheric 
wind  observations:  TID  horizontal  wavelength,  vertical 


wavelength,  intrinsic  frequency,  amplitude,  and  propagation 
direction  are  all  specified.  When  GPS-RI  receivers  are  co¬ 
located  with  other  instrumentation,  observed  and  calculated 
parameters  may  be  cross-validated.  Most  importantly,  TID 
horizontal  wavelengths  and  propagation  azimuths  may  be 
compared  to  GPS-RI  results  when  favorable  seeing  condi¬ 
tions  pertain,  providing  a  crucial  validation  of  this  technique. 
Calculated  values  of  intrinsic  frequency  and  vertical  wave¬ 
length  may  be  compared  with  ionosonde-derived  values. 

Low-latitude  TIDs  are  important  to  study  because  they 
are  a  key  mechanism  by  which  energy  and  momentum  are 
deposited  into  the  thermosphere.  Observation  of  TIDs  at 
altitudes  near  where  they  break  would  enable  a  direct  inves¬ 
tigation  of  this  process  for  the  first  time.  Additionally,  TIDs 
are  thought  to  be  a  key  seed  mechanism  of  plasma  depletions. 
While  evidence  of  TID  seeding  was  found  during  the  Spread  - 
FEx  campaigns  [37],  never  before  has  it  been  possible  to 
continuously  monitor  TID  activity  on  the  F-region  bottom- 
side  during  the  crucial  early-evening  hours,  when  seeding 
occurs.  An  investigation  of  the  presence/absence  of  TID 
activity  on  the  bottomside  will  help  determine  the  impor¬ 
tance  of  TIDs  relative  to  other  mechanisms  toward  seeding 
equatorial  plasma  depletions. 

Characterization  of  TID  large-scale  features  should  help 
establish  which  types  of  TIDs  are  most  effective  seeds.  For 
example,  is  a  minimum  perturbation  amplitude  needed  to 
seed  depletions?  TID  wavelength  and  propagation  azimuth, 
also  characterized  by  GPS-RI  arrays,  are  also  thought  to 
influence  seeding  effectiveness.  Observations  [38]  suggest 
that  TIDs  must  have  a  minimum  horizontal  wavelength  to 
effectively  seed  depletions.  Model  results  and  seeding  theory 
indicate  that  a  TID  must  propagate  nearly  parallel  to  the  geo¬ 
magnetic  equator  to  be  maximally  effective  at  seeding  plasma 
depletions.  GPS-RI  makes  it  possible  to  directly  investigate 
the  impact  of  the  variation  of  horizontal  wavelength  and 
propagation  azimuth  on  plasma  depletion  formation. 

7.  Conclusions 

We  have  presented  calculations  of  the  phase  velocity,  scale 
size,  and  angle  of  propagation  of  TIDs  that  are  likely  ori¬ 
ginated  by  gravity  waves  that  are  propagating  in  the  north¬ 
western  side  of  the  South  American  continent.  The  values 
that  were  obtained  using  the  SADM-GPS  and  the  CCM 
methods  are  consistent  with  values  of  medium-scale  TIDs 
observed  by  other  authors.  We  calculated  that  on  July  20, 
2008,  the  propagation  velocity  and  the  azimuth  direction  of 
the  TIDs  were  140  m/s  and  parallel  to  geographic  north. 

We  also  applied  the  CCM  method  to  stations  separated 
by  hundreds  of  km,  like  Huancayo,  Cuzco,  and  Piura.  It  was 
unexpected  to  find  a  large  coherence  between  the  TECP 
traces  from  all  the  adjacent  stations,  allowing  us  to  determine 
the  morphological  characteristics  of  the  TIDs.  This  result 
opens  the  possibility  of  using  measurements  of  the  TEC 
perturbations  by  the  GPS  receivers  already  deployed  in  South 
America  to  monitor  TID  activity.  We  have  conducted  two 
more  campaigns  using  a  small  network  of  three  receivers 
that  were  placed  near  a  new  ionosonde  that  operated  at 
the  Jicamarca  radar  site  between  2009  and  2010.  We  expect 
to  publish  these  results  in  the  near  future. 
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